The Quantum Emergence of Chaos 
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The dynamical status of isolated quantum systems, partly due to the linearity of the Schrodinger 
equation is unclear: Conventional measures fail to detect chaos in such systems. However, when 
quantum systems are subjected to observation - as all experimental systems must be - their dynamics 
is no longer linear and, in the appropriate limit(s), the evolution of expectation values, conditioned 
on the observations, closely approaches the behavior of classical trajectories. Here we show, by 
analyzing a specific example, that microscopic continuously observed quantum systems, even far 
from any classical limit, can have a positive Lyapunov exponent, and thus be truly chaotic. 
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There can be no chaos in the dynamics of isolated or 
closed quantum systems, a result which follows primar- 
ily from the linearity of the Schrodinger equation Q and 
the linear Hilbert space structure of the theory which, 
by virtue of the uncertainty principle, prevents the for- 
mation of fine-scale structure in phase space precluding 
chaos in the sense of classical trajectories. This leads 
to a widely recognized difficulty, as classical mechanics, 
which manifestly exhibits chaos, must emerge from quan- 
tum mechanics in an appropriate macroscopic limit 0. 
The key to the resolution of this apparent paradox lies 
in the fact that all experimentally accessible situations 
necessarily involve measured, open systems: the central 
importance of such situations in the context of chaos 
was first emphasized by Chirikov In a closely con- 
nected development, continuous quantum measurement 
theory Q has led to the successful understanding of 
the emergence of classical dynamics from the underly- 
ing quantum physics i i B i, and inequalities have 
been derived that encapsulate the regime under which 
classical motion, and thus classical chaos, exists 0. The 
transition to classical mechanics results from the local- 
ization of the quantum density matrix due to the infor- 
mation continuously provided by the measurement (it- 
self mediated by an environmental interaction), and the 
balancing of this against the unavoidable noise from the 
quantum backaction of the measurement. For a macro- 
scopic system, the Ehrenfest theorem holds as a result of 
localization and, simultaneously, the backaction noise is 
negligible, resulting in a smooth classical trajectory. 

While it has been established that observed quan- 
tum systems can be chaotic when they are macroscopic 
enough that classical dynamics has emerged, can they 
be chaotic outside this limit? This is the question we 
address here. By defining and computing the Lyapunov 
exponent for an observed quantum system deep in the 
quantum regime, we are able show that the system dy- 
namics is chaotic. Further, the Lyapunov exponent is not 



the same as that of the classical dynamics that emerges 
in the classical limit. Since the quantum system in the 
absence of measurement is not chaotic, this chaos must 
emerge as the strength of the measurement is increased, 
and we examine the nature of this emergence. 

The rigorous quantifier of chaos in a dynamical system 
is the maximal Lyapunov exponent The exponent 
yields the (asymptotic) rate of exponential divergence of 
two trajectories which start from neighboring points in 
phase space, in the limit in which they evolve to infinity, 
and the neighboring points are infinitesimally close. The 
maximal Lyapunov exponent characterizes the sensitivity 
of the system evolution to changes in the initial condi- 
tion: if the exponent is positive, then the system is expo- 
nentially sensitive to initial conditions, and is said to be 
chaotic. We apply this notion below to the observation- 
conditioned evolution of quantum expectation values. 

The evolution of a simple single-particle quantum sys- 
tem under an ideal continuous position measurement is 
given by the nonlinear stochastic master equation (SME) 
for the system density matrix (llj |: 



dp = — — [H, p] dt 



k[x, [x, p]]dt 
+^k{xp ^ px- 2{x)){dy - {x)dt) , 



(1) 



where the first term on the right hand side is due to uni- 
tary evolution, H being the Hamiltonian, and the sec- 
ond term represents diffusion from "quantum noise" due 
to the unavoidable quantum backaction of the measure- 
ment. The position operator is x, and the parameter k 
characterizes the rate at which the measurement extracts 
information about the observable, and which we will re- 
fer to as the strength of the measurement [l^ . The final 
term represents the change in the density matrix as a 
result of the information gained from the measurement. 
Here, dy is the infinitesimal change in the continuous 
output of the measuring device in the time dt. The con- 
tinuous output of the measuring device, y{t), referred 
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FIG. 1: Position distribution for the Duffing oscillator with 
measurement strengths k = 0.01 (red) and fc = 10 (green), 
demonstrating measurement-induced localization (fc = 10). 
The momentum distribution behaves similarly. 



to usually as the measurement record, is determined by 
dy = {x)dt + dW/^/%k where dW is the Wiener incre- 
ment, describing driving by Gaussian white noise |l3l| . 
The noise dW is due to the fact that the results of the 
measurement are necessarily random. (Note that the 
backaction and dW are uncorrelated with each other.) 
Thus on a given experimental run, the system will be 
driven by a given realization of the noise process dW . 
We will label the possible noise realizations by s. 

A single quantum mechanical particle is in principle 
an infinite dimensional system. However, for the pur- 
pose of defining an observationally relevant Lyapunov 
exponent, it is sufficient to use a single projected data 
stream: Here we choose the expectation value of the po- 
sition, {x(t)). The important quantity is thus the di- 
vergence, A(t) = \{x{t)) — (xgd(i))|, between a fiducial 
trajectory and a second trajectory infinitcsimally close 
to it. It is important to keep in mind that the system 
is driven by noise. Since we wish to examine the sensi- 
tivity of the system to changes in the initial conditions, 
and not to changes in the noise, we must hold the noise 
realization fixed when calculating the divergence. The 
Lyapunov exponent is thus 



A 



lim lim 

t^oo As(0)^0 



lnA,,(i) 



t 



lim \s{t) 



(2) 



where the subscript s denotes the noise realization. This 
definition is the obvious generalization of the conven- 
tional ODE definition to dynamical averages, where the 
noise is treated as a drive on the system. Indeed, un- 
der the conditions when (noisy) classical motion emerges, 
and thus when localization holds (Fig. it reduces to 
the conventional definition, and yields the correct clas- 
sical Lyapunov exponent. To combat slow convergence, 
we measure the Lyapunov exponent by averaging over an 
ensemble of finite-time exponents As {t) instead of taking 
the asymptotic long-time limit for a single trajectory. 



A key result now follows: In unobserved, i.e., isolated 
quantum dynamical systems, it is possible to prove, by 
employing unitarity and the Schwarz inequality, that A 
vanishes; the finite-time exponent, X{t), decays away as 
1/t This theorem codifies the expectation that, since 
the evolution is linear, any measure of chaos applied to it 
should yield a null result. As we have emphasized earlier, 
however, once measurement is included the evolution be- 
comes nonlinear and the Lyapunov exponent need not 
vanish. We now address this crucial question for a spe- 
cific example. 

The system we consider is the quantum Duffing os- 
cillator 'idl, which is a single particle in a double-well 
potential, with sinusoidal driving. The Hamiltonian for 
the Duffing oscillator is 



H = p^/2m + Bx^ - Ax^ + Ax cos(wt) 



(3) 



where p is the momentum operator, m the particle mass, 
and A, B and A determine the potential and the strength 
of the driving force. We fix the values of the parameters 
to be m = 1, B = 0.5, A = 10, A = 10 and uj = 6.07. 
The action of a system relative to h can be varied ei- 
ther by changing parameters in the Hamiltonian, or by 
introducing scaled variables so that the Hamiltonian re- 
mains fixed, but the effective value of h becomes a tun- 
able parameter. Here we employ the latter choice as it 
captures the notion of system size with a single number; 
the smaller U the larger the system size, and vice versa. 

To examine the emergence of chaos we will first choose 
fi — 10~^, which is small enough so that the system 
makes a transition to classical dynamics when the mea- 
surement is sufficiently strong. In this way, as we in- 
crease the measurement strength, we can examine the 
transformation from essentially isolated quantum evo- 
lution all the way to the (known) chaos of the clas- 
sical Duffing oscillator. To examine the emergence of 
chaos, we simulate the evolution of the system for k = 
5 X 10-^,10-3,0.01,0.1,1,10. When k < 0.01, the 
distribution is spread over the entire accessible region, 
and Ehrenfest's theorem is not satisfied. Conversely, 
for k = 10, the distribution is well-locahzed (Fig. 
and Ehrenfest's theorem holds throughout the evolution. 
Since the backaction noise, characterized by the momen- 
tum diffusion coefficient. Dp = h^k, remains small, at 
this value of k the motion is that of the classical system, 
to a very good approximation. 

Stroboscopic maps help reveal the global structural 
transformation in phase space in going from quantum to 
classical dynamics (Fig. The maps consist of points 
through which the system passes at time intervals sepa- 
rated by the period of the driving force. For very small 
k, (x) and (p) are largely confined to a region in the cen- 
ter of phase space. Somewhat remarkably, at fc = 0.01, 
although the system is largely delocalized, as shown in 
Fig. n nontrivial structure appears, with considerable 
time being spent in certain outer regions. By k — 1 the 
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FIG. 2: Phase space stroboscopic maps shown for 4 differ- 
ent measurement strengths, = 5 x fO~*, O.Of (top), and 
f, 10 (bottom). Contour lines are superimposed to provide 
a measure of local point density at relative density levels of 
0.05, 0.15, 0.25, 0.35, 0.45, and 0.55. 



localized regions have formed into narrower and sharper 
swirling coherent structures. At fc = 10 the swirls dis- 
appear, and we retrieve the uniform chaotic sea of the 
classical map (the small 'holes' are periodic islands). The 
swirls in fact correspond to the unstable manifolds of the 
classical motion. Classically, these manifolds are only 
visible at short times, as continual and repeated folding 
eventually washes out any structure in the midst of a 
uniform tangle. In the quantum regime, however, the 
weakness of the measurement, with its inability to crys- 
tallize the fine structure, has allowed them to survive: we 
emphasize that the maps result from long-time integra- 
tion, and are therefore essentially time- invariant. 

To calculate the Lyapunov exponent we implement 
a numerical version of the classical linearization tech- 
nique 15j, suitably generalized to quantum trajectories. 
The method was tested on a classical noisy system with 
comparison against results obtained from solving the ex- 
act equations for the Lyapunov exponents ,16j . The cal- 
culation is very numerically intensive, as it involves in- 
tegrating the stochastic Schrodinger equation equivalent 
to the SME over thousands of driving periods, and 
averaging over many noise realizations; parallel super- 
computers were invaluable for this task. 

We find that as t is increased, for nonzero fc, the value 
obtained for X{t) falls as 1/t, following the behavior ex- 
pected for fc = 0, until a point at which an asymptotic 
regime takes over, stabilizing at a finite value of the Lya- 
punov exponent as t oo. This behavior is shown in 
Fig.|3|for three different values of k. The Lyapunov expo- 
nent as a function of k is shown in Fig. 0] The exponent 




FIG. 3: Finite-time Lyapunov exponents X{t) for measure- 
ment strengths fc = 5 x 10~^, 0.01, 10, averaged over 32 
trajectories for each value of k (linear scale in time, top, and 
logarithmic scale, bottom; bands indicate the standard devi- 
ation over the 32 trajectories). The (analytic) 1/t fall-off at 
small k values, prior to the asymptotic regime, is evident in 
the bottom panel. The unit of time is the driving period. 



increases over two orders of magnitude in an approxi- 
mately power-law fashion as k is varied from 5 x 10~* to 
10, before settling to the classical value, Xci = 0.57. The 
results in Figs. 13 and 01 show clearly that chaos emerges 
in the observed quantum dynamics well before the limit 
of classical motion is obtained. 

We now compute the Lyapunov exponent for the quan- 
tum system when its action is sufficiently small that 
smooth classical dynamics cannot emerge, even for strong 
measurement. Taking a value of ?i = 16, we find that 
for fc = 5 X 10-^, A = 0.029 ± 0.008, for k = 0.01, 
A = 0.046 ± 0.01 and for k = 0.02, A = 0.077 ± 0.01. 
Thus the system is once again chaotic, and becomes more 
strongly chaotic the more strongly it is observed. From 
these results, it is clear that there exists a purely quantum 
regime in which an observed system, while behaving in a 
fashion quite distinct from its classical limit, nevertheless 
evolves chaotically with a finite Lyapunov exponent, also 
distinct from the classical value. 

It is worth pointing out that an analogous analysis can 
also be carried out for a continuously observed classical 
system. First one notes that an unobserved probabilistic 
classical system also has provably zero Lyapunov expo- 
nent: the average of x for an ensemble of classical par- 
ticles does not exhibit chaos, due to the linearity of the 
Liouville equation If we consider a noiseless ob- 

served chaotic classical system - possible since classical 
measurements are by definition passive (no backaction 



4 



' I Y~\ ity QED and nanomechanics experiments [ig. Experi- 

" mentally, one would use the known measurement record 

to integrate the SME (^); this provides the time evolution 
of the mean value of the position. From this fiducial tra- 
jectory, given the knowledge of the system Hamiltonian, 
the Lyapunov exponent can be obtained by following the 
procedure described here. 
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FIG. 4: The emergence of chaos: The Lyapunov exponent A 
as a function of measurement strength k. Error-bars follow 
those of Fig. 3, taken at the final time. 



noise) - then even the weakest meaningful measurement 
will, over time, localize the probability density, gener- 
ating an effective trajectory limit, and thus the classi- 
cal Lyapunov exponent, Ac; Noise can always be 
be injected into classical systems as an external drive, 
nevertheless, in the limit of weak noise, the system will 
once again possess the noiseless exponent Xci : In a clas- 
sical system the external noise is not connected to the 
strength of the measurement, so one can simultaneously 
have strong measurement and weak noise, which is pos- 
sible in the quantum theory only under specific condi- 
tions 0. As one way to understand this case, we can 
employ the quantum result as an intermediate step. Con- 
sider the quantum Lyapunov exponent at a fixed value 
of k (where A < Ac;) as in Fig. 01 If the value of h is now 
reduced, the dynamics of the system must tend to the 
classical limit as the quantum-classical correspondence 
inequalities of Ref. ^] are better satisfied. Thus the Lya- 
punov exponent in the classical limit of quantum theory 
- which, to a very good approximation, is just classical 
dynamics driven by weak noise - must tend to Ac;- If, 
however, the noise is not weak, an observed classical sys- 
tem, like a quantum system outside the classical regime, 
will also not be localized, and may well have an exponent 
different from Ac; . In addition, one may expect the non- 
localized quantum and classical evolutions to have quite 
different Lyapunov exponents, especially when h is large 
on the scale of the phase space, as quantum and classical 
evolutions generated by a given nonlinear Hamiltonian 
are essentially different The nature of the Lyaunov 
exponent for non-localized classical systems, and its re- 
lationship to the exponent for quantum systems is a very 
interesting open question. 

Finally, we emphasize that the chaos identified here is 
not merely a formal result - even deep in the quantum 
regime, the Lyapunov exponent can be obtained from 
measurements on a real system as in next-generation cav- 
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